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Abstract 

The flux of cosmic ray hadrons at tlie atmospheric depth of 820 g/ cm^ has been mea- 
sured by means of the EAS-TOP hadron calorimeter (Campo Imperatore, National 
Gran Sasso Laboratories, 2005 m a.s.l.). 

The hadron spectrum is well described by a single power law : 
Sh{Eh) = (2.25 ± 0.21 ± 0.34^2^'^) x 10-7(^)(-2-79±o.05) ^-'is~hi-^Ge\-^ 
over the energy range 30 GeV ^ 30 TeV. The procedure and the accuracy of the 
measurement are discussed. 



(*) Corresponding author. 
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The primary proton spectrum is derived from the data by using the COR- 
SIKA/QGSJET code to compute the local hadron flux as a function of the pri- 
mary proton spectrum and to calculate and subtract the heavy nuclei contribu- 
tion (basing on direct measurements). Over a wide energy range Eq = 0.5 -j- 50 
TeV its best fit is given by a single power law : 

S{Eo) = (9.8 ± 1.1 ± l.G'y) X 10-5(-g_){-2.80±o.o6) ^-2 ^-1 ^^-1 GeV-1. 

The validity of the CORSIKA/QGSJET code for such application has been 
checked using the EAS-TOP and KASCADE experimental data by reproduc- 
ing the ratio of the measured hadron fluxes at the two experimental depths 
(820 and 1030 g cm~^ respectively) at better than 10% in the considered 
energy range. 
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1 Introduction 

The spectrum of hadrons detected at different atmospheric altitudes retains 
significant information about the energy/nucleon spectrum of primary cosmic 
rays, which is dominated by the lightest component, i.e. the proton one. Its 
measurement has been carried on in the past, both at sea level [1-6] and at 
mountain altitude [7-11], using different experimental techniques, like emul- 
sion chambers, magnetic spectrometers and calorimeters. 
The knowledge of the primary proton spectrum is of main relevance for the 
understanding of the cosmic rays acceleration mechanisms and of the prop- 
agation processes in the Galaxy. Moreover, the proton component is mainly 
responsible for the uncorrelated particle production in the atmosphere: any 
uncertainties on the proton spectrum reffect in an uncertainty in the calcula- 
tion of the secondary particle ffuxes (tt and K) and thus, for example, on the 
knowledge of the atmospheric muon and neutrino ffuxes. A precise knowledge 
of such spectra is of particular importance to interpret the observational data 
from muon and neutrino detectors deep underground [12]. 
The measurement of the primary proton spectrum has been performed by 
means of satellite and balloon borne experiments [13-26] and indirectly de- 
rived by using ground based detectors [27-30]. In the region of tens of TeV, 
however, direct measurements lack statistics and moreover their energy deter- 
minations are not calorimetric and depend on the interaction parameters and 
their fluctuations. The data inferred from hadron measurements at ground 
level can therefore provide signiffcant new information. 

On the other side, the derivation of the information on the primary cosmic ray 
spectrum from hadron measurements, as well as the comparison of the results 
from different experiments, relies on the use of simulation tools describing 
the interaction and propagation of primary cosmic rays in the atmosphere. 
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The response of such hadron interaction models has therefore to be verified, 
especially considering that the recorded hadrons are the results of large fluc- 
tuations with respect to the average behavior. 

The EAS-TOP Extensive Air Shower array was located at Campo Impera- 
tore, 2005 m a.s.L, above the underground Gran Sasso Laboratories, with the 
aim of studying the cosmic ray spectrum in the energy range 10^^ -r- 10^^ eV 
through the detection of the different air shower components. 
In this paper, we present and discuss the results obtained in the study of the 
uncorrelated hadrons by means of the calorimeter of EAS-TOP, namely: 

a) the measurement of the hadron flux in the energy range 30 GeV 30 TeV; 

b) the derivation of the primary proton energy spectrum in the range 0.5 -j- 
50 TeV; 

c) the check of the propagation and interaction code (CORSIKA/QGSJET) 
used for the interpretation of the data. 



2 The detector and the hadron trigger 

The Muon and Hadron Detector of EAS-TOP is a 144 m^ calorimeter [31] 
made of nine layers, each composed by a 13 cm iron slab absorber (except for 
the uppermost plane which is unshielded), and three planes of limited streamer 
tubes, for a total depth of 818.5 g cm^"^. 

Two of the streamer tube layers (with 100 //m wire diameter at 4650 V voltage) 
act as tracking devices, and are read by a two-dimensional system based on the 
anode wires and external orthogonal pick up strips. The third one, which data 
are used for the present analysis, operates in saturated proportional mode (the 
wire diameter being 50 /xm, and the HV at 2900 V) for hadron calorimetry and 
EAS core studies. The signal charges are collected by a matrix of 840 (40x38) 
cm^ pads placed above the tubes; the pad signals are transferred to charge 
integrating ADCs with 15 bit dynamic range. The pad read-out is converted 
to the equivalent number of vertical particles by means of periodical calibration 
runs based on single muon triggers (pressure and temperature dependence of 
the induced charge being corrected for). 

Different sets of scintillators are placed in the apparatus for different aims; in 
particular, of the six ones lodged below the second absorber layer, three are 
used for hadronic trigger purposes. Each scintillator, of dimensions (80x80) 
cm^, is centered on a pad, viewed by two identical photomultipliers operating 
in coincidence and discriminated at the level of 30 m.i.p., corresponding to the 
energy loss of a 30 GeV proton incident on the calorimeter. They provide the 
"local hadron trigger" , which generates the read-out of the whole detector. 
For each scintillator a "tower" is defined, as the stack of 3 x 3 pads of the 8 
internal layers centered on the scintillator itself. The detector and its operation 
are fully described in ref. [31]; a scheme of a "tower" is shown in Fig.l. 



3 



I 2 Streamer tubes + 

' 1 Saturated PropDrtional 

J FeSlab 




Pad 

< === Tower === s 



Fig. 1. Schematic view showing the plane numbering and one of the defined "towers" 
inside the EAS-TOP calorimeter. 

3 Hadron selection, acceptance and energy calibration 

An event, recorded following the "local hadron trigger", is accepted as a 
hadron if: a) the cascade crosses at least three consecutive internal layers of 
the calorimeter, including the one positioned immediately below the trigger- 
ing scintillator, and b) the maximum energy release is recorded on the central 
pad of each plane of the corresponding "tower". This allows the selection of 
hadrons with energy above 30 GeV, and the definition of the angular accep- 
tance. 

The check of hadron selection, the detection efficiency, effective area, angular 
acceptance, and energy calibration have been obtained by means of simula- 
tions of the detector response based on the GEANT code [32] (with FLUKA 
option for hadronic cascades), including the full description of the apparatus. 
Protons at fixed primary energy and zenith angle have been generated and 
analyzed with the same procedure as the experimental data. 
Particular care has been put in the modelling of the chamber behavior in the 
saturated proportional mode; the saturation in the collected charge has been 
studied in detail and included in the simulation, as fully discussed in [31]. The 
modelling of the chamber response to large particle densities has been checked 
at the 50 GeV beam at CERN-PS, using a detector built by chambers with 
the same characteristics, read out and filling gas mixture as the ones oper- 
ating on site, but with length reduced to 3 m. Lead was used as absorber in 
front of the chambers in order to reach maximum particle density [31]. The 
chamber response was tested and found to agree with the model inside 2% 
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up to particle densities pch — 300 particles/cm^, corresponding to a 50 GeV 
electromagnetic shower after 4 cm of lead absorber. For iron absorber and 
the calorimeter geometry such particle density corresponds to hadrons with 
energy — 650 GeV. 

Such response, introduced in the quoted simulation, provides transition curves 
that can be directly compared with the experimental data. As shown in Fig.2, 
the difference between the two curves is always less than 10 % even at shower 
maximum, where the particle density is the highest, thus showing that the 
chamber behavior and saturation are well described at least up to 5 TeV (i.e. 
at particles densities at which the chamber response could not be directly 
tested) . 
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Fig. 2. Mean transition curves for hadrons in the calorimeter. Full circles: experi- 
mental data; empty circles: simulation. 



The verification of the hadron selection procedure has been performed by com- 
paring the shapes of the longitudinal developments for individual events with 
the expectations from the simulated ones (the agreement on the average tran- 
sition curves having been shown in Fig.2). For layers 1-7 (i.e. the ones shielded 
by more than two iron slabs, see Fig.l) the experimental and simulated Ni/Ntot 
distributions (i.e. the ratio between the equivalent particle number recorded in 
each layer and the total one in the tower) are in agreement inside the statistical 
errors (~ 10%). For layer 8 (shielded by a single iron slab), the contamination 
from the accompanying shower adds an excess of 15% of Ntot in 16% of the 
events, independent on A^^^. The effect does not alter the hadron selection 
and the spectrum measurement beyond the systematic effects discussed in the 
following. 

For the described triggering conditions, the effective area Aeff{Eh,9) was de- 
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termined using the same simulation code and taking into account the inef- 
ficiency of the trigger scintillators due to the 30 m.i.p. threshold. Such area 
includes the detection efficiency, which, concerning energy, rises above 65 % 
at ~ 130 GeV for vertical incidence inside the geometrical area of the cen- 
tral pad. As regards zenith angle, the efficiency at 30° is about 10 % of the 
vertical one. The selection condition therefore introduces a cut in the angular 
acceptance such that 90% of the events are found inside 22° from the vertical 
direction. 

The effective area of each "tower" is shown in Fig. 3 for 4 different zenith 
angles. 
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Fig. 3. Effective area Agff{Eh, 6) for each "tower" vs logioEf^ for 4 different hadron 
incidence angles. 

The hadron energy is inferred from the total charge induced on the 8 shielded 
layers of the defined "tower" (more than 95 % of the shower particles at all 
energies are contained inside a 20 cm radius from the hadron position). 
The conversion curve from the total number of particles induced in the "tower" 
(Ntot) to the primary hadron energy is shown in Fig.4. The energy resolution 
is a{Eh)/Efi = 15% at 1 TeV, worsening to 25% at 5 TeV due to leakage losses 
and to 30% at 30 GeV due to sampling losses. The dependence of the total 
number of particles on the hadron zenith angle is less than 3% up to 30°; the 
difference in the conversion curve between protons and pions impinging on the 
calorimeter is less than 2%. 

The possibility that the triggering and selection procedure includes more 
hadrons has been studied by means of a simulation of cascades in the at- 
mosphere through CORSIKA/QGSJET. It results that such hadron pile-up 
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effect, even at the highest energies {Eh >3 TeV), docs not alter the aver- 
age energy determination of more than 6%. As a test, to evaluate possible 
contaminations from the accompanying shower particles, the particle-energy 
conversion curve has also been obtained using the total charge induced on the 
five innermost planes only. No statistically significant difference was found in 
the hadron fiuxes obtained in the two cases over the considered energy range. 



4 The hadron flux 



About one million triggers were recorded in T= 14760 hours of effective live 
time used in the present analysis; 40832 of them survived the selection criteria 
and were classified as hadrons. 



The measured number of events in each energy bin for the fiux S{E, 9) is : 



il Eh+AEh 



rCr{Eu ^En + AE^) = j J S{E, e)TA,ff{E, e)dQdE (1) 

Eh 

The hadron fiux at zenithal angle 9 can be approximated as: 



S{E,e) = S{E)exp[ 



HE) 



(2) 
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Fig. 4. Total number of induced equivalent particles versus primary energy. 
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Table 1 



The measured hadron flux at 820 g/cn? . The given uncertainties are the statistical 
ones. 

where S(E) is the flux in vertical direction and x{6) is the atmospheric depth 
along 9. The attenuation length h-{E) has been derived using the CORSIKA 
[33] code to simulate the interactions and propagation of primary protons in 
air. In fact the hadron flux in the atmosphere includes both residual primaries 
and secondaries; at the EAS-TOP atmospheric depth, their ratio rises from 
~0.7 at 500 GeV to ~1.4 at 5 TeV. Therefore the obtained values of k{E) 
represent the full evolution of such mixture: A(£') ~ 114 g/cw? for QGSJET 
[34], and ~ 131 g/cm^ for HDPM [35], for the EAS-TOP altitude and range 
of zenith angles. 

Assuming a power law spectrum (7 = 2.7) inside each energy bin, the mean 
value Eh is obtained, the corresponding flux being S{Eh) = S{E){^)~'^ (a 
change of A7 = 0.1 in the spectral slope does not produce any appreciable 
difference in the resulting flux). 
The vertical flux is thus : 



gr^ ^ ^ rCr{En ^Eh + AE,) 

27rTEl J J E-^exp[ - ^]A,ff(E, e)sinededE 

The recorded number of events and the experimental hadron fluxes at the 
atmospheric depth of 820 g/cvo? are listed in Tab.l with the corresponding 
statistical uncertainties. 
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The following sources of systematic uncertainties have to be considered: 



1. The uncertainty in the evaluation of the effective area, — 12% at all 
energies. 

2. The uncertainty in the hadron angular distribution, which reflects in the 
evaluation of the attenuation length A.{E). A comparison between two dif- 
ferent models (QGSJET and HDPM) in the CORSIKA frame shows that 
the differences in A(E) reflect in a flux uncertainty ^ ~ 5%. 

3. The uncertainty in the energy assignment to each single hadron, for a spec- 
tral slope 7 ~ 2.7, results in a flux uncertainty ^ ~ 7%. This value reaches 
10% at the highest energies, as shown by the comparison of the measured 
and simulated longitudinal shower proflles. 

4. An uncertainty ^ ~ 15% on the flux, due to the different behavior and 
efficiency of the triggering scintillators and to the different calibrations and 
stability of the corresponding "towers". 

A total systematic energy dependent uncertainty ^ ~ 15% is obtained from 
the first 3 items. To this, the 15% constant systematic uncertainty due to item 
4 has to be added. 

The hadron flux is fltted by a power law from 30 GeV up to 30 TeV as 

ShiEn) = (2.25 ± 0.21) x 10-^(^)(-2-^9±o.o5) m'^s-'sr-'GeV-' (4) 
with = 0.91 and is shown in Fig.5. 

In the fltting procedure (and in the plot), the energy dependent systematic 
uncertainties have been included; the 15% energy independent systematic ef- 
fect has to be added. 

The hadron flux is compatible, within the errors, with a single power law. This 
has been tested by performing the same flt in independent narrower energy 
ranges, the resulting slopes being shown in Fig.6. 



5 The primary proton spectrum 

The primary proton spectrum is derived from the data by: 

a) checking the hadron propagation code in the atmosphere; 

b) subtracting from the measured hadron spectrum the contribution of heavy 
primaries; 

c) minimizing the difference between the measured and the expected hadron 
fluxes on the basis of different primary proton spectra. 

a) The region of interest coincides with the energy range in which QGSJET 
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Fig. 5. The hadron flux at 820 g/cm?. The best fit (4) is shown superimposed to the 
data. 
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Fig. 6. Slopes of the independent fits to the hadron flux. The fitting energy range is 
shown. The circle and dashed line show the slope as found in (4). 

(the hadronic interaction model used to describe the cosmic ray interaction and 
propagation in the atmosphere) has been directly checked against accelerator 
data [37,38], both concerning the leading particle and the secondary produc- 
tion physics. Its reliability to reproduce the present data has been checked by 
comparing its predictions to the measured ratio of hadron fluxes at sea level 
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E,(GeV) 

Fig. 7. Ratio of the experimentally measured hadron fluxes by EAS-TOP and KAS- 
CADE (squares) compared to the expectation if the proton+helium primary spectra 
by J A CEE ( circles ) or R UN JOB ( stars ) are assumed. 

(KASCADE [36], 1030 g/cm^) and mountain altitude (EAS-TOP, 820 g/cm^). 
Primary protons and helium nuclei were generated in quasi vertical direction 
6 < 5°, with energy spectra according to JACEE and RUN JOB [25,26] and 
the expected hadron fluxes at each observation level were calculated. As shown 
in Fig. 7, the expected ratio does not depend on the differences between such 
primary spectra, and it is compatible with the measured one within the sta- 
tistical uncertainties, the comparison leading to a x^=l-2/d.o.f. On average, 
the model reproduces the experimental ratio at better than 10%. 
We remind that the general features of the model relevant for the calculation 
of the hadron flux (and therefore object of the test) are the combination of the 
total cross section and inelasticity concerning the contribution of the surviving 
primaries, and the very forward production cross section for the contribution 
of the secondaries. We therefore conclude that QGSJET, as implemented in 
CORSIKA, can be reliably applied in the considered energy range in the de- 
scription of the uncorrelated hadron fluxes at different atmospheric depths 
and therefore can be applied between the top of the atmosphere and the EAS- 
TOP observation level, thus allowing to derive the primary nucleon flux from 
the present measurement. 

b) The contribution to the hadron flux from helium primaries has been evalu- 
ated using their spectrum as directly measured by the balloon experiments. In 
order to derive the systematic uncertainties of the procedure, both the RUN- 
JOB {'jHe = 2.80) and JACEE {jHe = 2.68) data have been used and their 
contribution subtracted from the experimental hadron flux. At Eh — 1 TeV, 
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such contribution is 15% and 29% from RUNJOB and JACEE respectively; 
the heavier nuclei one is less than 10 %. 

c)The primary proton spectrum is obtained as the one minimizing the differ- 
ence between the measured hadron spectrum (after subtraction of the Helium 
contribution by means of the afore described procedure) and the expected one 
from simulated proton primaries. Extensive simulations have been carried on, 
generating primary protons in quasi vertical direction {6 < 5°), with energy 
extracted on power law spectra with slope varying between 2.5 and 3.2. The 
number of simulated events is such that the number of hadrons in each en- 
ergy bin be much higher than the experimentally collected one. Most of the 
contribution to each hadron energy bin comes from different primary energy 
regions; hadrons of energies in different ranges, e.g. = 0.1 -r 0.2, 0.2 
0.5, 0.5 -r- 1, 1 ^ 2, 2 5, > 5 TeV are produced by primaries with median 
energy Emed — 0.5, 2, 4, 10, 20, 55 TeV respectively. The data thus allow to 
get information on the primary proton spectrum in the range 0.5 ^ 50 TeV. 
Assuming a primary spectrum of the power law form S{Eq) = SqEq^ , the 
normalization factor 5*0 and the slope 7 have been obtained minimizing the 
differences between the calculated and the measured number of hadrons in 
each energy bin. The minimizations have been carried on by taking into ac- 
count both the statistical and the energy dependent systematic uncertainties 
in the hadron flux. 

The data are well described by power law spectra in the energy range 0.5 
-j- 50 TeV, with best fits, for the case of subtraction of the RUNJOB Helium 
spectrum: 

S{Eo) = (1.05 ±0.16) X 10-4(^)(-2-80±0.05) ^-2 g-lgj.-l Q^y-l 

and for the case of subtraction of the JACEE Helium spectrum: 
S{Eo) = (0.91 ±0.15) X 10-4(^)(-2-80±o-06) ni-^ s-'sr-' GeV"! 

Including the 7% uncertainty in the helium contribution and the 15% constant 
systematic uncertainty on the measured hadron flux into a global systematic 
error term, the result can be summarized as follows : 

S{Eo) = (9.8 ± 1.1 ± X 10-5(^)(-2-80±o.o6) m-^s-isr-iQeV-i (5) 

The obtained proton spectrum is shown in Figs. 8 and 9. The full area and 
the shaded lines (in the two figures respectively) include the systematic and 
statistical uncertainties of the measurement. 
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Fig. 8. Primary proton spectrum; the full area represents the result of this mea- 
surement and includes the systematic and statistical errors. Results from different 
experiments are shown for comparison. The straight line represents a fit from [40J. 
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Fig. 9. Primary proton spectrum multiplied by Eq'^ . Symbols as in Fig. 8. The hori- 
zontal line represents the present result, Eq.5, with errors. 
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6 Conclusions 



The hadron flux has been measured over a wide energy range (30 GeV^30 
TeV) by means of the EAS-TOP hadron calorimeter at the atmospheric depth 
of 820 g/ cm^. The spectrum is weU described by a single power law in the whole 

range : 

S{Eh) = (2.25 ± 0.21 ±0.34^2^^) x IQ-^ {^)(-^-^^±om) ^-2 3-1 ^^-i Q^y-i^ 
Taking into account the contamination from heavier nuclei, on the basis of 
direct measurements, the primary proton spectrum is obtained between 0.5 
and 50 TeV and is found to be compatible with a single slope power law: 
S{Eo) = (9.8 ± 1.1 ± X 10-5(^)(-2-80±o.o6) ^^^-2 g-i g^-i GcV-^. 

A systematic uncertainty of about 7 % due to the uncertainty in the helium 
flux is included. The data match very well with the direct measurements over 
a wide energy range, usually not available to a single experiment, where direct 

measurements become statistically poor. 

The reliability of the CORSIKA/QGSJET interaction and propagation code, 
which is used to propagate the hadrons in the atmosphere and to compute the 
heavy nuclei contribution, is directly checked in this energy range by compari- 
son with accelerator data and, concerning the direct application to the present 
measurement, through its capability to reproduce the ratio of hadron fluxes as 
measured at two different atmospheric depths by EAS-TOP and KASCADE, 
at 820 and 1030 g/cm^ respectively. 
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